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We outline an algorithm for an efficient solution of the Quantum Counting problem using Adi- 
abatic Quantum Computation (AQC). We show that using local adiabatic evolution, a process in 
which the adiabatic procedure is performed at a variable rate, the problem is solved as efficiently 
as the analogous circuit-based algorithm, i.e., quadratically faster than the corresponding classical 
algorithm. The above algorithm provides further evidence for the potentially powerful capabilities 
of AQC as a paradigm for efficient problem solving on a quantum computer. We also discuss the 
algorithm in relation to the problem of Quantum Phase Detection. 
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I. INTRODUCTION 

The paradigm of Adiabatic Quantum Computation 
(AQC), proposed by Farhi et al. [l[ about a decade ago 
is a simple yet intriguing approach to problem solving 
on a quantum computer. Unlike the leading paradigm 
of circuit-based quantum computing, AQC is an ana- 
log, time-continuous method which does not require the 
design and use of quantum gates. As such, it can in 
many ways be thought of as a simpler and perhaps more 
profound method for performing quantum computations 
while also being easier to design and implement experi- 
mentally @,[|. 

The potential advantages of AQC over the usual 
circuit-based and other computation paradigms have nat- 
urally generated substantial research, and intensive ef- 
forts have been made towards finding computational 
problems that are solved more efficiently on an adiabatic 
quantum computer (henceforth, a quantum annealer) 
than on a classical computer (see, e.g., Refs. [3-ll2j and 
references therein). One of the most important results 
in that context, which serves as an additional indication 
to the potential capabilities of the AQC paradigm, is the 
existence of a theorem stating that adiabatic quantum 
computing is polynom ially equivalent to circuit model 
quantum computation p~3l [T^| . 

Despite the above equivalence, it is not yet understood 
in which cases adiabatic quantum computation is just as 
efficient as circuit model quantum computation (or even 
potentially polynomially more efficient). It has therefore 
become of great interest to look at problems for which 
circuit model algorithms are known, and to see whether 
equally-efficient adiabatic algorithms can be found. To 
date, however, there are almost no examples to show that 
AQC is as efficient as the circuit-model paradigm. One 
notable exception is the adiabatic Grover's search algo- 
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rithm devised by Roland and Cerf [Tq . demonstrating 
that the quadratic speed-up over classical algorithms first 
shown by Grover |l6[ using circuit-based computation, 
can also be gained with adiabatic quantum algorithms 
(we discuss the algorithm in more detail later). 

In what follows, we consider another problem for 
which, much like Grover's search problem, there is a 
circuit-model algorithm that is quadratically faster than 
its classical counterpart. This is the problem of 'Quan- 
tum Counting' (T3, IU| in which one has to estimate the 
number of occurrences of a specific item in an unsorted 
database. We show that the Quantum Counting prob- 
lem may also be solved on a quantum annealer in an 
efficient manner via the use of Local Adiabatic Evolu- 
tion (LAE) [ill]. We shall see that the efficiency of the 
proposed method is on par with the analogous circuit- 
based algorithm and so is quadratically faster than the 
corresponding classical algorithm. 

The paper is organized as follows. In the next section 
we briefly discuss the principles of the Quantum Adia- 
batic Algorithm that is the heart of AQC. We then de- 
scribe is some detail the concept of Local Adiabatic Evo- 
lution (LAE) in Sec. IIIII In Sec. IIV1 we analyze the dy- 
namical equations of the QAA for the adiabatic Grover's 
algorithm. These are later used in Sections [V] and IVII to 
devise an efficient quantum-adiabatic algorithm for solv- 
ing the problem of Quantum Counting. In Sec. IVIII we 
briefly discuss the suggested algorithm in the context of 
quantum phase detection, and finally, in Sec. I Villi we 
conclude with a few comments. 



II. QUANTUM ADIABATIC ALGORITHM 
(QAA) 

We first describe the Quantum Adiabatic Algorithm 
(QAA) which provides the general approach for solving 
problems by casting them in the form of optimization 
problems 

Within the framework of the QAA approach, the solu- 
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tion to an optimization problem is encoded in the ground 
state of a Hamiltonian H p . To find the solution, the 
QAA prescribes the following course of action. As a 
first step, the system is prepared in the ground state 
of another Hamiltonian Hd, commonly referred to as the 
driver Hamiltonian. The driver Hamiltonian is chosen 
such that it does not commute with the problem Hamilto- 
nian and has a ground state that is fairly easy to prepare. 
As a next step, the Hamiltonian of the system is slowly 
modified from to H pi using the linear interpolation, 



i.e. 



H{s) =sH p + (l- s)H d , 



(1) 



where s(t) is a parameter varying smoothly with time, 
from at t ~ to 1 at the end of the algorithm, at 
t = T. If this process is done slowly enough, the adia- 
batic theorem of Quantum Mechanics (see, e.g., Refs. (l9j 
and [20[) ensures that the system will stay close to the 
ground state of the instantaneous Hamiltonian through- 
out the evolution, so that the system finally arrives at 
a state that is close to the ground state of H p . At this 
point, measuring the state will give the solution of the 
original problem with high probability. 

In the simple case where s(t) varies from zero to one 
at a constant rate, the runtime T must be chosen to be 
large enough so that the adiabatic approximation holds: 
this condition determines the efficiency, or complexity, of 
the QAA. A condition on T can be given in terms of the 
instantaneous eigenstates ||m)} and eigenvalues {E m } of 
the Hamiltonian H(s), as fa. |21| 



r > 



1 max s Vbi(g) 
e min s g 2 (s) 



where e is a small number, inversely proportional to the 
running time of the algorithm, that controls the expected 
probability for failure of the algorithm, g(s) is the first 
excitation gap Ei(s) — E (s) and Vbi(s) = |(0|dff/ds|l)| 
(in our units h = 1). 



III. LOCAL ADIABATIC EVOLUTION 

As is well established by now, LAE can provide a signif- 
icant speed-up in calculation times if one utilizes the con- 
cept of 'Local Adiabatic Evolution' (LAE), in the frame- 
work of which the adiabatic parameter s(t) is varied not 
at a constant rate but rather at a variable rate, slowing 
down in the vicinity of the minimum gap and speeding- 
up in places where the gap is large . The procedure is 
based on formulating a 'local' Landau-Zcncr condition for 
each value of the adiabatic parameter s. This condition 
can be written as: 



ds 



< € 



Voi(s) ' 



(3) 



where g(s), Vbi(s) and e are as in Eq. @. 



One example in which one may utilize the concept 
of LAE and which in turn demonstrates the potentially 
powerful capabilities of Adiabatic Quantum Computa- 
tion was first illustrated by Roland and Cerf jl5l |22| . 
In this example, an efficient adiabatic algorithm for solv- 
ing Grover's search problem [l6| which describes a search 
for a marked item in an unstructured database, is found. 
More specifically, Roland and Cerf showed that while the 
application of the adiabatic theorem with a linear s(t) 
results in a running time that is of order N - N being 
the number of items in the database - a carefully cho- 
sen variable rate of s(t) yields a running time that scales 
like y/N, i.e., a quadratic speed-up is gained, similarly 
to the original result by Grover, that was found for the 
circuit-based model. 

LAE however can only be efficiently used in cases 
where one has proper knowledge of the exact behavior of 
the gap and relevant matrix elements of the system for 
the problem at hand. This is normally not the case. The 
amenability of Grover's search problem to LAE stems 
from the fact that the gap of the system (as well as the 
relevant matrix element) can be calculated explicitly be- 
forehand, even if the specific ground state is not known. 
With a proper choice of problem and driver Hamiltonians 
(which are discussed next), the gap in the Grover search 
problem is given by the simple expression: 



g(.m s) = 4(1-^(1-*) , 



(4) 



where rj = M/N is the ratio of the number of solutions 
M to the size of the database N. Moreover, the matrix 



element 



(dH/ds 



can also be calculated in a straight- 



(2) forward manner to give: 



Vbi(s) 



dH/ds 



VvO- - rj) 



(5) 



The condition for local adiabatic evolution for the 
Grover problem is thus simply given by: 



< e- 



(6) 



A less strict inequality using the fact that Vqi(s) < 1 
yields the more convenient bound of (l5j : 



(7) 



Integrating in time either of the two equations above 
between t = (also, s = 0) and t = T (s=l) yields 
immediately: 



1 



Tex -y/N/M 



1 



(8) 



in the limit of A — > oo, i.e., a quadratic speed-up over 
any classical search algorithm, a result also shared by 
Grover's original circuit- based algorithm [l6j . 
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IV. ANALYSIS OF THE ADIABATIC GROVER 
PROBLEM 

An important property of the gap of Grover's search 
problem, Eq. (j4j, may be formulated as follows. 

9{Vl\ s) < g(V2; s) for any 771 < 772 , (9) 

where 771 and 772 are two arbitrary ratios. The above prop- 
erty leads to the interesting (and somewhat surprising) 
attractive feature of LAE for the Grover search problem: 
If one runs the algorithm as if the ratio of solutions is 77* 
even though the actual ratio of solutions 77 may not be 
known, then one still obtains at the end of the adiabatic 
run a state that is very close the symmetric superposition 
of all solution states, as long as 77* < 77. This is because 
choosing the adiabatic profile with 77*, namely, 

^ = eff 2 (77,; S ), (10) 

also satisfies the inequality Eq. ([7]) for any 77* < 77. 

The above feature is to be contrasted with the circuit- 
based Grover's search algorithm, in which case one must 
know in advance the number of solutions in order to ob- 
tain a solution to the problem with a very high probabil- 
ity. Within the circuit model, the problem of not know- 
ing the number of solutions is resolved by a preliminary 
application of a different algorithm known as Quantum 
Counting designed to estimate the number of solutions 
M prior to the application of Grover's search algorithm. 

In what follows we show that LAE, with the aid of 
the unique properties of the adiabatic Grover's problem 
presented above, may be used to solve adiabatically the 
problem of Quantum Counting, i.e., the problem of esti- 
mating the number M of occurrences of a specific item 
within an unsorted database of size N. 

Before doing so however, we first analyze in some de- 
tail the dynamics of the adiabatic algorithm for Grover's 
search problem (a similar derivation is briefly discussed 
in Ref. [HI]). Specifically, we derive here the equations 
describing an unsorted-database search with an unknown 
number (cquivalcntly, ratio) of solutions M (or 77) using 
an adiabatic profile s(t) with a ratio 77*. We then cal- 
culate the probability of ending up in a solution state in 
two limiting cases. 

A. Dynamics of the Grover-adiabatic problem 

Consider the problem Hamiltonian for the Grover 
problem [l5j : 

H p = l- \m)(m\. (11) 

Here the M states meM are eigenstates of the compu- 
tational basis with eigenvalue zero (all other states have 
eigenvalue one). The problem Hamiltonian is therefore a 



diagonal matrix of size N with M zero entries and N — M 
entries of one along the diagonal. The driver Hamiltonian 
is simply given by: 

it d = -\4>){4>\, (12) 

where \<p) is the fully-symmetric state: 

v i— 1 

Obviously, the ground-state of the driver Hamiltonian, 
i.e., the state that the system is prepared in for the adi- 
abatic evolution, is 

m = o)> = \4>) ■ (14) 

Applying the Schrddinger equation 

i±m))=&w)) as) 

to = (c (rj),Ci(t),...,Ci(rj)),...,Cjv(i)) simplifies 

greatly due to symmetry considerations: The state of 
the system has only two distinct amplitudes. The ampli- 
tudes Ci(t) that correspond to 'solution' states and those 
that correspond to 'non-solution' states. Let us denote 
the former set by a(t) and that of the latter by b(t). 
The Schrddinger equation thus becomes the two coupled 
equations: 

ia = -(1 - s) [770+ (1 - 77)6] (16) 
ib = -(1 - s) [rja + (1 - 77)6] +sb, (17) 

Note that N, the size of the problem and M, the number 
of solutions appear in the equations only as the ratio 
77 = M/N. 

The time-dependence of a and b is easily transformed 
into s dependence (as there is a one-to-one correspon- 
dence between the two in our case) . The time-derivative 
transforms in the usual way: d/dt = (ds/dt) x d/ds, and 
the two coupled equations now become: 

ieg 2 (r)*;s)a' = -(1 - s) [rja + (1 - rj)b] (18) 
ieg 2 (rj*;s)b' = -(1 - s) [rja + (1 - 77)6] +sb . (19) 

Here, the prime symbol (') stands for differentiation with 
respect to the new independent variable s and we have 
used Eq. (|TU)) for ds/dt where the ratio 77* may in general 
be different than 77. As already noted earlier, a choice 
of 77* < 77 results in a true adiabatic evolution and is 
expected to yield at the end of the adiabatic process, a 
state that is very close to the equal superpositions of all 
solution states. 

In addition, the adiabatic process dictates the initial 
conditions 

a{s = 0) = b(s = 0) = l, (20) 

which correspond to the fully-symmetric state (i.e., an 
equal superposition of all possible classical states). The 
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specific value of 1 chosen here leads to the simple nor- 
malization condition 

r,\a\ 2 + (l- V )\b\ 2 = l, (21) 

that, similarly to the equations of motion, does not ex- 
plicitly depend on N or M but only on their ratios. 



B. End probabilities 

At the end of the evolution, at s = 1, the final state of 
the system is a linear combination of the equal superpo- 
sition of solutions states, namely, of oc X)me.M \ m ) 
and the equal superposition of the non-solution states 
iVWsoi) oc YlmiM \ m )- The more the evolution is adi- 
abatic, the larger the amplitude of the solution-states 
superposition will be. While the exact values of the two 
complex components at s = 1 can not be computed an- 
alytically, one can obtain an asymptotic solution in two 
important limits. One is the limit of small e, i.e., when 
the allowed error in the final state goes to zero. In the 
e — >• limit, the two equations above fall under the WKB 
category |23| , and have a series solution of the form 

oo 

a(s) = e^W withcr(s) = ^(ie) fe a & (s) (22) 

fc=0 

which can directly be plugged into the equations [note 
that b(s) is a simple linear combination of a(s) and <z'(s)]. 
Solving the resulting equations for the growing powers of 
e term by term, one may calculate the probability to 
obtain a solution state, i.e., the squared overlap between 
the final state at s = 1 and the optimal solution |"0soi): 

Psoi = \{i>(s = l)\i> sol }\ 2 (23) 

= hm 'Wl . 

.->i»?|o( S )P + (l-,7)|&(*)| a 

Calculation of the above series solution leads to the fol- 
lowing asymptotic solution for P so \ in the limit of small 
e: 

P soI (e^0) = l-e 2 r){l - v ) + 0(e 4 ). (24) 

To leading order in e, the probability to find a solution 
does not depend on value of 77* , as long as one chooses a 
small enough value for e. As expected, when e goes to 
(i.e., when no error is allowed), the probability to end up 
in the equal superposition of solution states approaches 
1 (in this case, the running time of the algorithm goes to 
infinity) . 

The more interesting limit however, is the one in which 
e is small but kept fixed while N, the dimension of the 
Hilbert space (alternatively, the size of the database) goes 
to infinity. This limit translates to the limit where the 
ratios rj = M/N and rj* — M*/N go to zero. Here, we as- 
sume that the number of solutions M docs not scale with 



N . This can be imposed by setting an iV-indcpcndent up- 
per bound on the number of solutions, M max . (We note 
that in the case where the number of solutions M scales 
with N, the Quantum Counting problem becomes less 
interesting, as the classical probability of finding a solu- 
tion, pjciassicai) _ ^yy^y _ ^ becomes of order unity). 

In the limit of N — > 00, following a careful asymp- 
totic analysis, we obtain for P so i, quite expectedly, the 
following 'Landau-Zener type' probability: 

TV 77 7T M 

P so i(N ^ 00) w 1 -e~^ = 1 -e"*^ . (25) 

As expected, the probability of obtaining a solution in- 
creases with decreasing e and also, as discussed in the 
previous section, with decreasing M*, which is set by 
the adiabatic profile function. Most importantly, in the 
N — > 00 limit, the above probability has a well estab- 
lished ./V-indcpcndcnt value. The above probability will 
become useful in what follows when we address the prob- 
lem of Quantum Counting. 

V. QUANTUM COUNTING: ZERO VS 
NONZERO NUMBER OF SOLUTIONS 

Before addressing the full-complexity of the general 
quantum counting problem, let us first address the spe- 
cial case of finding whether or not the number of solutions 
M (number of zero entries in the problem Hamiltonian) is 
zero. Circuit-based algorithms giv e a running time that 
is on the order of 0(VN) [13, Ell! that is, a quadratic 
speedup over the analogous classical algorithm. 

A successful adiabatic algorithm may be constructed 
in the following manner. One runs the QAA algorithm 
twice: Once with the original problem Hamiltonian of 
Hp, Eq. (fTTj) . and and a second time with H' p = 1 — H p . 
The adiabatic profile function in both cases is chosen to 
be such that s(t) obeys ds/di = eg 2 (rj = 1/N,s), i.e., 
suitable for efficiently solving the case of any nonzero 
number of solutions. 

Now, if the number of solutions M is zero, then nei- 
ther of the problem Hamiltonians will have any effect on 
the state of the system, and the state will remain the 
initial state \4>) throughout the adiabatic evolution, ex- 
cept for an undetectable global phase. A measurement of 
the z-magnetization at the end of the run will produce a 
random configuration with energy (for H p ) in one run 
and energy 1 (for H') in the other. 

If on the other hand the number of solutions is nonzero, 
the Grover LAE procedure ensures us that at the end of 
each of the two adiabatic runs, at s = 1, the final state of 
the system will be very close to \ip so i), the equal super- 
position of zero-energy states with a probability greater 
than P so \ = 1 — e~T^ (in the N — > 00 limit), and a mea- 
surement along the z direction will produce a zero-energy 
state with this probability in both runs. 

We conclude then, that two runs of the adiabatic al- 
gorithm whose running time is T oc \/N are sufficient 
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for establishing whether or not the number of solutions 
is zero. 



VI. QUANTUM COUNTING: THE GENERAL 
CASE 



The general case of estimating M can be formulated 
based on the probability P so i given in Eq. (J2SJ). Let us 
first briefly analyze the classical analogue of that same 
problem. In the classical case, a random sampling of 
items in the database will give p(°j assloal > = M/N = r\ for 
each picked item. The number of successes after k runs 
will therefore be distributed binomially with the unbiased 
mean of ^( classlcal ) — 77 an d a sampling error of: 



At) 



(classical) _ V Vi 1 ~ V) _ VV 

Vk ~ Vk 



(26) 



In terms of the number of solutions M , we find that the 
error in M is 



^J^J (classical) 



MN 



(27) 



i.e., for the error to scale as y/M, the number of trials 
k required for the search, has to scale linearly with the 
dimension of the database N. 

Going back to the corresponding quantum problem, 
the question becomes: How many independent runs k of 
the adiabatic procedure, each run producing a solution- 
state with probability P so i [as given by Eq. (pj)) ]. are nec- 
essary to yield an unbiased estimate of M with an error 
rate that, similarly to the classical case, scales like \fM~1 
Noting that the probability generated by each adiabatic 
run, does not depend on the size of the database ./V (in 
contrast to the probability in the classical quick 
calculation analogous to the one done for the classical 
case above yields: 



AM 



4eM* 



1 g4eM* 



1 



(28) 



The number of trials needed to yield an error that scales 
like y/~KI is therefore: 



16e 2 M? e 



k oc 



- 1 



Mtt 2 



(29) 



Choosing to be somewhat larger than M max (for ex- 
ample, 10 times the cutoff value) , the exponent above can 
be approximated by taking only the two leading terms of 
a Taylor-series expansion, and we end up with the re- 
quired number of trials k oc M* . 

Note that AM, while it may depend on M, does not 
scale with N. The number of trials needed to obtain 
any desired error is therefore N- independent. Now, the 
running time of each trial on the other hand scales like 



y/N/M* [see Eq. ©]. The overall running time of the 
quantum algorithm therefore scales like y/~N (as M* docs 
not scale with N) and the algorithm is therefore quadrat- 
ically faster that its classical analogue. 

VII. PHASE DETECTION 

The Quantum Counting algorithm outlined above, 
may also be used to solve phase estimation problems 
within which one is required to estimate the eigenphases 
(i.e., the phases of the eigenvalues) of a given unitary 
matrix. The above quantum counting algorithm may be 
utilized here in the sense that the (unnormalized) fully- 
symmetric initial wave function \<j)) may be re-expressed 
in the form: 



|0) = y/rj\a) + y/T^vi 
where we have defined: 

\a) = - V \m) , 
77 £ — ' 

1—77 



(30) 

(31) 
(32) 



Replacing 77 = M/N with the appropriate angle of 
sin 2 (8/2) reveals that estimating M is the same as es- 
timating the angle 9 for the wave function: 



sin (9/2) \a) + cos (9/2) |/3) , 



(33) 



meaning that quantum counting may be thought of as 
a phase estimation algorithm, where the phase here is 
represented by the angle 6. 



VIII. SUMMARY AND CONCLUSIONS 

We have given a prescription for solving the Quantum 
Counting problem in which the number of occurrences 
M of a specific item within an unsorted database of size 
N is to be estimated. We have shown that using Local 
Adiabatic Evolution (LAE) on Grover-type problems, the 
running time of the algorithm scales like VN, i.e., it is 
quadratically faster than the corresponding classical al- 
gorithm and on par with the circuit-based result. 

It would be interesting to see whether the algorithm 
outlined above can be used as the basis for the solution 
to the more general problem of Quantum Phase Esti- 
mation of which the phase detection problem discussed 
previously may be thought as a special case, and within 
which one is required to estimate the eigenphasc of an 
eigenvector of a given unitary operator. 

The importance of the above results, aside from the 
implications of having an efficient quantum-adiabatic al- 
gorithm for this and related problems, lies in the fact that 
it is one of very few examples that illustrate the power 
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and potential that are encompassed in adiabatic quantum 
computation, providing an adiabatic quantum algorithm 
that is as efficient as its circuit-model counterpart. 

Acknowledgments 

We thank Peter Young, Eleanor Rieffel and Jeremie 
Roland for useful comments and discussions. We ac- 



knowledge partial support by the National Security 
Agency (NSA) under Army Research Office (ARO) con- 
tract number W911NF-09-1-0391, and in part by the 
National Science Foundation under Grant No. DMR- 
0906366. 



[1] E. Farhi, J. Goldstone, S. Gutmann, J. Lapan, A. Lund- 
gren, and D. Preda, Science 292, 472 (2001), a longer ver- 
sion of the paper appeared in arXiv:quant-ph/0104129. 

[2] L. M. K. Vandersypen, M. Steffen, G. Breyta, C. S. Yan- 
noni, M. H. Sherwood, and I. L. Chuang, Nature 414, 
883887 (2001). 

[3] Z. Bian, F. Chudak, W. G. Macready, L. Clark, and 

F. Gaitan (2012), (arXiv:1201.1842). 
[4] T. Hogg, Phys. Rev. A 67, 022314 (2003). 
[5] E. Farhi, J. Goldstone, and S. Gutmann (2002), 

(arXiv:quant-ph/0201031). 
[6] E. Farhi, J. Goldstone, S. Gutmann, and D. Nagaj, Inter- 
national Journal of Quantum Information 6, 503 (2008), 

(arXiv:quant-ph/0512159). 
[7] B. Altshuler, H. Krovi, and J. Roland (2009), 

(arXiv:0908.2782v2). 
[8] S. Knysh and V. N. Smelyanskiy (2011), 

(arXiv:1005.3011). 
[9] A. P. Young, S. Knysh, and V. N. Smelyanskiy, Phys. 

Rev. Lett. 101, 170503 (2008), (arXiv:0803.3971). 
[10] A. P. Young, S. Knysh, and V. N. Smelyanskiy, Phys. 

Rev. Lett. 104, 020502 (2010), (arXiv:0910.1378). 
[11] I. Hen and A. P. Young, Phys. Rev. E. 84, 061152 (2011), 

arXiv:1109.6872v2. 



[12] I. Hen, Phys. Rev. E. 85, 036705 (2012), 

arXiv:1112.2269v2. 
[13] D. Aharonov, W. van Dam, J. Kempe, Z. Landau, 

S. Lloyd, and O. Regev, SIAM Journal of Computing 

37, 166 (2007). 

[14] A. Mizel, D. A. Lidar, and M. Mitchell, Phys. Rev. Lett. 
99, 070502 (2007). 

[15] J. Roland and N. J. Cerf, Phys. Rev. A 65, 042308 (2002). 

[16] L. K. Grover, Phys. Rev. Lett. 79, 325 (1997). 

[17] G. Brassard, P. Hoyer, and A. Tapp, Automata, Lan- 
guages and Programming 1443, 820 (1998). 

[18] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, England, 2000). 

[19] T. Kato, J. Phys. Soc. Jap. 5, 435 (1951). 

[20] A. Messiah, Quantum Mechanics (North-Holland, Ams- 
terdam, 1962). 

[21] G. H. Wannier, Physics 1, 251 (1965). 

[22] J. Roland and N. J. Cerf, Phys. Rev. A 68, 062312 (2003). 

[23] C. M. Bender and S. A. Orszag, Advanced Mathematical 
Methods for Scientists and Engineers (McGraw-Hill, New 
York, USA, 1978). 



